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We review some difficulties of the standard picture of confining fields which are 
viewed usually as 'mild', or quasiclassical field configurations. The paradoxes are 
cn : naturally resolved by recent lattice observations of selftuned fluctuations which 

exhibit both ultraviolet and infrared scales. The ultraviolet scale is provided by 
CO ■ the lattice spacing a and is manifested in the non-Abelian action density associated 

with the fluctuations while the physical scale is exhibited by their geometrical 
pH , characteristics. The data might suggest existence of a dual description of Yang- 

p») ■ Mills theories on the fundamental level. 



1. Introduction 



Vacuum fields are responsible for the confinement. This dogma is rarely put 
in doubt. However, much less is known specifically on the corresponding 
field configurations. Still the common belief is that the confining fields are 
non-perturbative a and 'mild', see, e.g., a review 2 and references therein. 
The evidence is , however, indirect. Moreover, in the Minkowski space it is 
difficult actually to formulate what 'soft' field means because the vacuum 
should look the same in all the inertial frames. 

The problem seems to be much more tractable in the Euclidean space. 
One is inclined to think of a single size characterizing one or another fluc- 
tuation. Moreover, there are two intrinsic scales in QCD, ultraviolet and 
infrared. We will have in mind the lattice regularization and by the ultra- 
violet cut off will understand, therefore, the lattice spacing a. The infrared 
scale is provided by Aqcd- 

There are well known examples which demonstrate relevance of these 
scales to the vacuum fluctuations. First, consider zero-point fluctuations. 
They have no intrinsic scale built in and the typical size is determined by 



a See, however, 



1 
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an external probe. For example, the vacuum expectation value of the gluon 
condensate (that is, a point-like probe) is given by: 

<o|(g^) 2 io> « {N 'j 1} , a) 

where G a is the gluonic field strength tensor and N c is the number of 
colors, and is clearly ultraviolet-dominated. The very fact of the ultraviolet 
divergence in the matrix element (1) is not specific, of course, for the lattice 
and well known from calculations of the vacuum energy density in any 
formalism. The lattice allows to unambiguously fix (1). 

On the other hand, the standard image for non-perturbative fluctuations 
is provided by instantons. The typical size of the instantons is of order 
(Aqcz?) -1 - Moreover there exist detailed models of the instanton liquid 
developed along these lines 3 . 

Thus, the standard picture is that at short distances zero-point fluctu- 
ations dominate while the confinement is due to soft quasiclassical fields, 
something like instantons. This picture can be probed, for example, through 
the sum-rule technique 4 and its generalizations and claims many phe- 
nomenological successes, see, e.g., 3 and references therein. 

However, there has been accumulating evidence that this standard pic- 
ture misses sometimes even large, leading effects 5 . We will briefly review 
the issue in Section 3. In Section 4 we argue that the problems of the stan- 
dard picture are resolved if one takes into account the newly discovered 
branes 6 . The branes are two-dimensional surfaces whose total area scales 
in the physical units while the action density is of perturbative order (1). 
The branes represent a new kind of selftuned fluctuations which unify the 
scales a and Aqcd- This, written version of the lectures presents only an 
overview of the actual content of the lectures given at the Summer Institute. 
Further details and references can be found in other talks of the author 7 . 



2. Difficulties of the standard picture 

Confinement 

We are considering pure Yang-Mills theory and the criterion of con- 
finement is the existence of a linear potential at large distances for heavy 
external quarks. The lattice data can indeed be fitted by a very simple 
form, 

const 
~R 

at all the distances. 



V QQ( R ) « "-5- + a - R > ( 2 ) 
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The point crucial for our discussion here is that the instanton liquid 
model does not reproduce the linear piece at all 8 . Thus, the successes 
of the model in other cases turn into a problem for the standard picture. 
Indeed, the model seems to be the best realization of the standard picture 
but still misses the confinement. 

Current correlators 

A standard way of probing vacuum fields is provided by measuring cur- 
rent correlators, < 0\j(x), i(0)|0 > where j(x) are local currents constructed 
from the quark and gluon fields, see, e.g., 4 ' 3 . The strategy is to start from 
short distances, x <C Aqc D and approach 'intermediate' distances where 
non-perturbative effects become sizable. In the most simplified form the 
theoretical predictions look as: 

<0|j(x),j(0)|0> « (parton model) (l +c,x 4 - < 0|(G^) 2 |0 > so/t +...) , 

(3) 

where the coefficient Cj depends on the current considered and calcula- 
ble while < 0|(G^„) 2 |0 > S oft stands for the contribution of the soft non- 
perturbative fields, < 0|(G£„) 2 |0 > so /t~ ^qcd- I R actual applications, 
there can be other vacuum condensates representing the non-perturbative 
fields. Another variation of (3) is an explicit calculation of a single- instanton 
contribution 3 . Also, one commonly keeps the first-order perturbative con- 
tribution of order a s (x) as well. 

In many cases, the approximations like (3) work well. However, there 
are cases when (3) fails, see in particular 5 . Moreover, there is accumulat- 
ing evidence that in cases when the Born (parton-model) approximation is 
dominated by gluon propagator the leading correction is of order 

<0|j(a;),j(0)|0> « (parton model) (1 + b j A 2 QCD -x 2 + ...) , (4) 

for the latest example of this kind and references see 9 . 

There is no model-independent way to relate the coefficients bj in (4) 
in various channels since the quadratic correction is not captured by the 
operator product expansion. Phenomenologically, however, the model with 
a non- vanishing short-distance gluon mass m g (see second paper in Rcf 5 ) 
turns to be successful. Within this model one modifies the propagator of a 
gluon by replacing 

1/q 2 - 1/q 2 + \m\\l ( ? , (5) 

where q is the momentum of the gluon, q 2 ^> m 2 . It is worth emphasizing 
that the model (5) applies in the Born approximation when higher order 
corrections in a s are ignored. (Modification of the Born approximation 



February 1, 2008 22:49 WSPC/Trim Size: 9in x 6in for Proceedings 



SI03 



4 

valid at all q 2 has been introduced in 10 ). In particular, the Cornell potential 
(2) at short distances exhibits a correction induced by (5). Indeed, the Born 
approximation in this case is given by one-gluon exchange and this is just 
the case when one expects the dominance of the quadratic correction. 

3. Long perturbative series 
3.1. Expectations 

The difference between 'long' and 'truncated' perturbative series is crucial 
in view of the theorem 11 that terms of order x 2 ■ Aq CD are calculable 
perturbatively. Let us explain this point in more detail. 

Generic perturbative expansion for a matrix element of a local operator 
looks as: 

oo 

( O } — (parton model) ■ (l + a n a™ ) , (6) 

n=l 

where we normalized the anomalous dimension of the operator O to zero. 
Note also that a s is the bare coupling, a s « 1. 

In fact, expansions (6) arc only formal since the coefficients a n grow 
factorially at large n: 

\a n \ ~ c?-n! , (7) 

where Cj are constants. Moreover, there are a few sources of the growth (7) 
and, respectively, Cj can take on various values, for review see, e.g., 12 . The 
factorial growth of a n implies that the expansion (6) is asymptotic at best. 
Which means, in turn, that (6) cannot approximate a physical quantity to 
accuracy better than 

A ~ exp ( - 1/aa.) ~ (xQCD-a) > ( 8 ) 

where bo is the first coefficient in the /^-function. To compensate for these 
intrinsic uncertainties one modifies the original expansion (6) by adding the 
corresponding power corrections with unknown coefficients. 

In case of the gluon condensate the theoretical expectations can be 
summarized as: 

(0 | ^i)( G ^) 2 |0) « a s (N \ 1] a n a: + [const^ ■ k% c D ) , 

n—l 

(9) 
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where 



b a s 



and terms proportional to J^q CD correspond to < 0|(G^„) 2 |0 > so ft, see (3). 
It is the latter quantity which enters the QCD sum rules, see 4 . 

A conspicuous feature of the prediction (9) is the absence of a quadratic 
correction, compare (4). However, Eqs (4) and (9) are not necessarily in 
contradiction with each other. The point is that the phenomenological 
analysis relies on short or even very short truncated perturbative series 
plus a power-like correction while the theorem (9) keeps a long perturbative 
series (plus the first power-like correction) 



3.2. Numerical results 

Numerically, this problem was studied in greatest detail in case of the gluon 
condensate on the lattice 13 . In terms of the lattice formulation the gluon 
condensate is nothing else but the average plaquette action. The result can 
be summarized in the following way. Represent the plaquette action (P) 
as: 

(P) « P£ rt + b N a 2 K% CD +c N a 4 A 4 QCD , (10) 

where the average plaquette action (P) is measurable directly on the lat- 
tice and is known to high accuracy, Pp ert is the perturbative contribution 
calculated up to order N: 

Ppert = 1 - £W" , (11) 
n=l 

and, finally coefficients 6/v,cjv are fitting parameters whose value depends 
on the number of loops N. Moreover, the form of the fitting function 
(10) is rather suggested by the data than imposed because of theoretical 
considerations. 

The conclusion is that up to ten loops, N = 10 it is the quadratic 
correction which is seen on the plots while c N are consistent with zero. 
However, the value of b N decreases monotonically with growing N 13 . The 
factorial divergence (7) is not seen yet and perturbative series reproduces 
the measured plaquette action at the level of 10~ 3 . Finally, at the level 
10~ 4 the Aq CD term seems to emerge 13 . 
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4. Branes 

4.1. Coexistence of two scales 

We have already mentioned that instantons do not confine. A natural 
question is then, what are the confining configurations. Actually the answer 
is also more or less known. Namely it is commonly believed (for a recent 
review and references see 14 ) that it is the monopolcs and central vortices 
that confine. The next question is then, why we did not account for these 
fluctuations, say in (9) or (3). The answer is pure technical, at first sight. 
Namely, the monopoles and vortices are defined in terms of projected, not 
original non-Abelian fields. The projection, in turn, is determined non- 
locally and the track to the original Yang-Mills fields is lost. However, the 
general expectation seems to be that the monopoles and vortices correspond 
to soft, bulky fields and as far as dependence on Aqcd is concerned are not 
much distinguishable from the instantons. 

These expectations turned to be not true. Namely both the monopoles, 
see 15 and references therein, and vortices 6 appear to be associated with 
an excess of the non-Abelian action which is divergent in the ultraviolet: 

(S mon ) ~ (S vort ) « 0.54-^ , (12) 

where L is the length of the monopole trajectory, A is the area of the vortex. 

It is most remarkable that the infrared scale is also relevant to the 
branes. Namely, it has been known since some time (for review see 14 ) 
that the densities of monopoles and vortices scale in physical units. The 
corresponding densities are defined as: 

Lperc — ^PpercVi , A vor t — GpvortVi , (13) 

where V4 is the volume of the lattice and L pers is the total length of the 
percolating cluster while A vor t is the area of the vortices. The densities 
(13) scale in physical units and are independent on the lattice spacing. 
According to the latest measurements: 

Pperc = 7.70(8) fm- 3 , Kort « 4.0(2) /m" 2 , (14) 

see 16 and 6 , respectively. Moreover, the monopole trajectories lie on the 
P-vortices 17 ' 6 . 

4.2. Selftuning 

Naively, one would expect that monopoles and vortices with action (12) 
propagate only very short distances, L ~ a, A ~ a 2 . However, both 
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monopoles and vortices form clusters which percolate through the whole 
of the lattice volume V4. This implies cancellation of the ultraviolet diver- 
gences between action and entropy. This cancellation is easy to quantify 
in case of monopoles (particles). Namely the propagating mass, m of the 
monopole is in fact not the radiative mass M(a), where M(a) = S mon /L, 
see (12). The relation between the two masses is as follows: 



const / . . In 7 



m 



2 

raon if 



( M (a) - — ) . (15) 



The In 7 term here is due to the entropy (see, e.g., 18 ). Observation (14) 
implies cancellation between the action and entropy. Moreover, there is no 
parameter to tune to ensure this cancellation. Thus, the monopoles are 
rather selftuned. The same is true for the vortices. Both the tension and 
entropy are ultraviolet divergent 6 . 



4.3. Quadratic correction to the gluon condensate 

Monopoles and vortices are defined, for a given configuration of the vac- 
uum fields, for the whole of the lattice. Thus, they are seen as a nonlocal 
structure. Moreover, they are manifestly non-perturbative. Indeed, the 
probability 6(plaq) for a particular plaquette to belong to a brane has been 
found to be proportional to: 

9(plaq) w (const) cxp( - l/b g 2 (a)) ~ (a ■ A QC d) 2 ■ (16) 

On the other hand, the branes have an ultraviolet divergent tension which 
assumes a kind of locality. 

To make contact with the continuum theory it is useful to evaluate con- 
tribution of the branes into local or quasi-local matrix elements. The gluon 
condensate turns to be the easiest case. Indeed, combining Eqs (14) and 
(12) one gets for the contribution of the vortices to the gluon condensate: 

( (G^) 2 ) V ort » 0.3 GcV 2 a- 2 , (17) 

We see that the mysterious a 2 ■ Aq CD correction, see Sect. 3.2., gets its 
explanation in terms of the branes. The mixture of the the two scales, a 
and Aqcd, exhibited by this correction appears to be a manifestation of 
selftuning of the branes. 
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5. Status of theory 

5.1. Why branes? 

The lattice data on the monopoles and vortices have been accumulating 
since long. However, they were analyzed almost exclusively within the con- 
finement problem. Discovery of the ultraviolet divergences, see (12), makes 
the challenge to the theory much more direct. Indeed in an asymptotically 
free theory one should be able to understand short distances from first 
principles. At present, however, the phenomenology is by far ahead of the 
theory. 

Observation of the branes is gratifying from the theoretical point of view. 
Indeed, the results (12) imply that at least at presently available lattices the 
size of the monopoles is not resolved. Moreover, gluons are already (approx- 
imately) free particles at such distances. On the other hand, it is well known 
that there exists no consistent field theory with both electric (color) and 
magnetic dynamical charges. Existence of the branes implies that the mag- 
netic and electric (color) charges are separated in space. The color charges 
live in the bulk while the magnetic charges live on the two-dimensional 
branes (which percolate through the four- dimensional Euclidean space). 

5.2. Constraints from asymptotic freedom 

The asymptotic freedom implies that at short distances the only degrees of 
freedom relevant are those of free gluons. How do we count the degrees of 
freedom? Usually through ultraviolet divergences. The best known example 
is the /3-function which counts logarithmic divergences. The logarithmic 
divergences are singled out since they are independent on details of the 
cut off. However, on the lattice the power divergences are also uniquely 
defined. And they are much easier to study. In particular, if we think in 
terms a scalar complex field 4>m describing monopoles then we can introduce 
a vacuum expectation value {\4>m\ 2 ) b - For an elementary monopoles we 
would have (|^m| 2 ) ~ o, 2 - This is not allowed since we may not have new 
particles at short distances. What is allowed is 

(I0m| 2 ) ~ A 2 QC D , (18) 
and this is a constraint from the asymptotic freedom. 



b The condensate can be thought of as non-perturbative part of the gauge invariant 
condensate of dimention two, see 19 . 
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Remarkably enough, this constraint can be rewritten 7 in terms of the 
total length of the monopole trajectories (in (13) we considered only the 
percolating cluster while now arc including finite clusters as well c ): 

p^ on < (const)^^ . (19) 

This constraint turns to be satisfied by the data 16 . On the other hand, 
from pure geometrical point of view (19) implies that the monopoles ac- 
tually percolate on a two-dimensional surface, not on the whole of the 
four-dimensional space. Which means branes. 

In other words the data on the monopole action, see (12) plus asymptotic 
freedom of the original YM theory imply existence of branes. 

5.3. Manifestations of the duality 

As we mentioned above, the original formulation of the YM theory does 
not explain the selftuning of the monopole action, compare (12) and (15). 
Indeed the geometrical factor In 7 depends on the type of the lattice (cubic) 
and is meaningless in the continuum. Imagine, however, that there were a 
theorem proving existence of the dual formulation in the continuum limit. 
Then the selftuning would be derived since it is a necessary condition for 
the monopoles to survive in the continuum limit. 

Also, compare the contribution to the gluon condensate from instantons 
and monopoles. The instantons are added to the perturbative expansion, 
see (9). The contribution of the branes, to the contrary, cannot be added 
to the perturbation theory but is, instead, dual to the perturbative series, 
see Sects. 3.2, 4.3. The reason is that the branes 'belong' to the dual world. 

Acknowledgments 

I am thankful to the organizers of the Summer Institute, and especially to 
Prof J. Kubo and Prof H. Terao, for the invitation and hospitality. The 
paper was worked out while the author was visiting the Yukawa Institute 
for Theoretical Physics of the Kyoto University. I am thankful to Prof. T. 
Kugo and other members of the group for the hospitality. 

References 

1. V.N. Gribov, Eur. Phys. J. CIO, 71 (1999). 



c For further applications of the percolation theory see , in particular, 2() - 



February 1, 2008 22:49 WSPC/Trim Size: 9in x 6in for Proceedings 



SI03 



2. Yu. Dokshitzor, "QCD phenomenology", hep-ph/0306287. 

3. V. Shuryak, Th. Schafer, Ann. Rev. Nucl. Part. Set. 47, 359 (1997). 

4. M.A. Shifman, A.I. Vainshtein, V.I. Zakharov, Nucl. Phys. B147, 385 (1978). 

5. V.A. Novikov, M.A. Shifman, A.I. Vainshtein, V.I. Zakharov, Nucl. Phys. 
B191, 301 (1981); 

K.G. Chetyrkin, S. Narison, V.I. Zakharov, Nucl. Phys. B550, 353 (1999), 
(hep-ph/9811275); 

F.V. Gubarev, M.I. Polikarpov , V.I. Zakharov, "Physics of the power correc- 
tions in QCD", hep-ph/9908292; 

E. V. Shuryak, "Probing the boundary of the nonperturbative QCD by small 
size instantons", hep-ph/9909458. 

6. F.V. Gubarev, A.V. Kovalenko, M.I. Polikarpov, S.N. Syritsyn, and V.I. Za- 
kharov, Phys. Lett. B574, 136 (2003) , hep-lat/0309034; 

A.V. Kovalenko, M.I. Polikarpov, S.N. Syritsyn and V.I. Zakharov, "Interplay 
of monopoles and P-vortices", hep-lat/0309032. 

7. V.I. Zakharov, "Hidden mass hierarchy in QCD", hep-ph/0204040; 
" 'Branes' in the lattice SU(2) gluodynamics", hep-ph/0306361; 

"Fine tuning in lattice SU(2) gluodynamics vs continuum-theory constraints", 
hep-ph/0306362; 

"Nonperturbative match of the ultraviolet renormalon" , hep-ph/0309178; 
"Hints on dual variables from the lattice SU(2) gluodynamics", hep- 
ph/0309301. 

8. D. Chen, R.C. Brower, J.W. Negele, E.V. Shuryak, Nucl. Phys. Proc. Suppl. 
73, 512 (1999), (hep-lat/9809091). 

9. S. Narison, V.I. Zakharov, Phys. Lett. B522, 266 (2001). 

10. G. 't Hooft, "Perturbative confinement", hep-th/0207179. 

11. M. Beneke, V. I. Zakharov, Phys. Rev. Lett. 69, 247 (1992). 

12. M. Beneke, Phys. Rept. 317, 1 (1999); 

R. Akhoury, V.I. Zakharov, Nucl. Phys. Proc. Suppl. 54A, 217 (1997) (hep- 
ph/9610492). 

13. G. Burgio, F. Di Renzo, G. Marchesini, E. Onofri, Phys. Lett. B422 219 
(1998), (hep-ph/9706209); 

F. Di Renzo and L. Scorzato, JHEP 0110, 038 (2001), (hep-lat/0011067) 

R. Horsley, P.E.L. Rakow, and G. Schierholz, Nucl. Phys. Proc. Suppl. 106, 
870 (2002), hcp-lat/01 10210. 

14. J. Greensite, Prog. Part. Nucl. Phys. 51, 1 (2003) , hep-lat/0301023. 

15. V.G. Bornyakov, et al, Phys. Lett. B537 291 (2002). 

16. V.G. Bornyakov, P.Yu. Boyko, M.I. Polikarpov, V.I. Zakharov, Nucl. Phys. 
B672, 222 (2003). 

17. J. Ambjorn, J. Giedt, J. Greensite, JHEP 0002, 033 (2000). 

18. J. Ambjorn, "Quantization of geometry" , hep-th/9411179. 

19. F.V. Gubarev, L. Stodolsky, V.I. Zakharov, Phys. Rev. Lett. 86, 2220 (2001). 

20. M.N. Chernodub, V.I. Zakharov, Nucl. Phys. B669, 233 (2003) (hep- 
th/0211267); 

M.N. Chernodub, K. Ishiguro, K.Kobayashi, T.Suzuki, "Numerical determi- 
nation of monopole entropy in pure SU(2) QCD", hep-lat/0306001. 



